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We work out the basic analysis on dynamics near QCD critical point (CP) by dynamic renormal- 
ization group (RG). In addition to the RG analysis by coarse graining, we construct the nonlinear 
Langevin equation as a basic equation for the critical dynamics. Our construction is based on the 
generalized Langevin theory and the relativistic hydrodynamics. Applying the dynamic RG to the 
constructed equation, we derive the RG equation for the transport coefficients and analyze their 
critical behaviors. We find that the resulting RG equation turns out to be the same as that for the 
liquid-gas CP except for an insignificant constant. Therefore, the bulk viscosity and the thermal 
conductivity strongly diverge at the QCD CP. We also show that the thermal and viscous diffusion 
modes exhibit critical slowing down with the dynamic critical exponents zthormai ~ 3 and Zviscous ~ 2 
, respectively. In contrast, the sound propagating mode shows critical speeding up with the negative 
exponent Zsound ~ —0.8. 

PACS numbers: 12.38.-t, 05.70.Jk, 25.75.Nq 



00 



Oh. 
I 

CD 



> 
in 
oo 

o 



X 



I. INTRODUCTION 

An interesting feature of the phase diagram in quan- 
tum chromodynamics (QCD) is the possible existence of 
the critical point (CP), which is predicted by various ef- 
fective models of QCD and suggested by lattice QCD 
simulations. The critical point is the end point of the 
first order phase transition line existing in the low tem- 
perature (T) region Then, the significance of this 
QCD CP is that the phase transition at this point is of 
second order, and thereby we can expect critical phenom- 
ena due to large fluctuations of various quantities at this 
point 

Then, a fundamental problem arises; what is the soft 
modes of the QCD CP? A hint is that the baryon-number 
susceptibility Q diverges at the CP as first suggested 
in and subsequently demonstrated in It has 

been established that the fluctuating modes of conserved 
densities are the soft modes at the CP [H, ■ Although 
the a mode seems to be the soft mode, it couples to 
the density fluctuation at flnite density ^] and remains 
massive 

Furthermore, as a critical phenomenon, some authors 
suggested a divergent behavior of bulk viscosity at the 
QCD CP lisl, although their validity of their argument is 
very much controversial for instance, the ansatz 

for the spectral function adopted in [Tcj may not neces- 
sarily be true [TTj], and a microscopic calculation by the 



'Electronic address: y-minami@ruby.scphys.kyoto-u.ac.jpl 
^ In fact, the QCD matter seems to have extremely rich structure 
in the phase diagram with one or multiple critical points ^jj- 
and even accompanied with inhomogeneous phases at low T 
01, although the very existence of a CP may be questioned 
according to [3, 0| . 

Such a fast mode as the a near the QCD CP is called a slaved 
mode Il5|l , because its slow dynamics is controlled by the density 
fluctuation for the QCD CP. 



relativistic Boltzmann equation [19| shows that the bulk 
viscosity is finite at the CP. Thus, it is still not obvi- 
ous whether the transport coefficients will diverge at the 
QCD CP. 

In fact, as is known in condensed matter physics, the 
critical divergence of the transport coefficients is a com- 
mon phenomenon at a CP, such as at the liquid-gas CP, 
and originate from a universal mechanism; nonlinear fluc- 
tuations of macroscopic variables cause the divergence 
[l^H^I. This implies that microscopic processes, as may 
be described by such as the Boltzmann equation, would 
give only a minor contribution to the critical divergence 
of these quantities, if an y. T he dynamic renormalization 
group (RG) theory 0,1111 is a standard technique for 
critical dynamics, which systematically incorporate the 
macroscopic fluctuations causing the divergent behavior 
of transport coefficients. In this theory. We must con- 
struct a nonlinear Langevin equation as a basic equation 
for the critical dynamics. The construction goes as fol- 
lows. First, the slow variables are identified for describing 
the critical dynamics. Next, the thermodynamic poten- 
tial for the slow variables is constructed to determine 
the static property of the system. Finally, the stream- 
ing terms, causing the dynamic-nonlinear coupling, and 
the kinetic coefficients are determined for respectively de- 
scribing time reversible and irreversible changes of the 
slow variables. We note that the streaming term is ab- 
sent in the simple Brownian motion. 

The general theory of the critical dynamics as de- 
scribed above tells us that an essential ingredient is to 
properly construct the nonlinear Langevin equation for 
the critical dynamics. As far as we know, this is the 
first attempt for the QCD CP. Our construction of the 
Langevin equation is based on the generalized Langevin 
theory, by Mori (22I . [29| , and the relativistic hydro- 
dynamics, because the slow variables are identified as 
lon g- wa velength fluctuations of the conserved densities 
013 HSl ; we construct the streaming terms from con- 
tinuity equations and the potential condition, which is a 
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general condition for streaming terms (20|, |2l|. Also, we 
use the thermodynamic potential for the 3d Ising system 
as that for the QCD CP because the static universality 
class is the same as 3d Ising class [H, [Hill. Finally, 
we determine the kinetic coefhcients from a relativistic 
hydrodynamic equation, here the so called Landau equa- 
tion |24| taken. In consequence, we shall show that the 
Langevin equation differs from it for the liquid-gas CP 
due to relativistic effects, although the dynamic univer- 
sality class of the QCD CP is conjectured as of the liquid- 
gas CP [Hill. 

After such construction, we apply the dynamic RG to 
the Langevin equation and derive the RG equations for 
the transport coefficients. Consequently, to our surprise, 
these RG equations turn out to be the same as for the 
liquid-gas CP except for a irrelevant constant, although 
the Langevin equations are different. Therefore, the bulk 
viscosity and the thermal conductivity strongly diverge 
and can be more important than the shear viscosity near 
the QCD CP. We shall also show that the thermal and 
viscous diffusion modes exhibits critical slowing down, 
whereas the sound mode critical speeding up. 

This paper is organized as follows. In Sec. II, we shortly 
review the theory of critical dynamics. In Sec. Ill, we con- 
struct the nonlinear Langevin equation for the QCD CP 
by the generalized Langevin theory and the relativistic 
hydrodynamics. In Sec. IV, we analyze the critical ex- 
ponents of the transport coefficients and dynamic criti- 
cal exponents by the dynamic RG. The final section is 
devoted to a summary and concluding remarks. In Ap- 
pendixes El |B] and we give the detailed derivation of 
the RG equations for the transport coefficients. 



II. REVIEW OF CRITICAL DYNAMICS 



Since the present work is based on the theory of critical 
dynamics in condensed matter physics (20l . [2l| , we now 
shortly review it for self-containedness. 



A. Critical divergences of transport coefficients 

The critical divergence of transport coefficients (or dif- 
fusion constants) is a common phenomenon, for instance, 
at the critical point of the liquid-gas, ferromagnetic tran- 
sitions and so on [2^ [2l| . The important point is that 
the critical divergence originates from a universal mech- 
anism; nonlinear fluctuations of macroscopic variables 
causes the divergence [13, [13] ■ 

Here, we briefly illustrate how macroscopic nonlinear 
fluctuations cause the critical divergence, taking the ther- 
mal conductivity near the liquid-gas CP as an example 
[28| . The thermal conductivity is given by the Kubo for- 
mula as follows. 



where ^(r, t) and T are the heat current and temperature, 
respectively, and (• • •) denotes the statistical average in 
the thermal equilibrium state. The heat current g(r, t) is 
supplied from two sources: one is due to a microscopic 
process as calculated by a microscopic theory, such as the 
Boltzmann equation, and the other is due to nonlinear 
fluctuations of macroscopic variables j29| : 



ro + 



(2) 



where (/micro and gmacro respectively denote the micro- 
scopic and macroscopic currents. The macroscopic pro- 
cess causing the heat current is identified as the entropy 
density convected by fluid velocity fluctuation. Thus, we 
have 

^macro ^ 

6sSv, (3) 

where Ss and dv respectively denote the fluctuations of 
the entropy density and the fluid velocity. The macro- 
scopic current Eq. ^ is of the second order in fluctua- 
tions and hence negligible far from the CP. However, it 
becomes the dominant part near the CP, since the fluc- 
tuations are enhanced there. We see that Eq. ([l} now 
has the following form 



A = An 



+ / dr 



dr 



dt{q{r,t)q{0,0)), 



(1) 



dt{Ss{r, t)Sv{r, t)Ss{0, 0)Sv{0, 0)), 

(4) 

where Amicro is the thermal conductivity coming from the 
microscopic current. Recalling that the entropy density 
fluctuation is a soft mode near the liquid-gas CP, we see 
that the second term of Eq. (|4]) diverges at the CP. This 
is the mechanism causing the critical divergence of trans- 
port coefficients. 

Let us call the transport coefficients, such as Amicro, 
coming from microscopic processes the bare transport 
coefficients, and those including the contributions from 
the nonlinear macroscopic fluctuations the renormalized 
ones. Then, we need not to study the critical divergence 
of transport coefficients by a microscopic theory since 
the divergence originates from only the macroscopic pro- 
cesses. The dynamic RG [H [U, [2!, [3l[, [12] is the stan- 
dard theory incorporating such nonlinear macroscopic 
fluctuations. In this theory, we must construct a nonlin- 
ear Langevin equation as a basic equation for the critical 
dynamics. 



B. Generalized Langevin equation 

We first note that if the dynamic variables are divided 
into slow and fast ones, the slow dynamics can be well 
described by a Langevin equation. We stress that such 
the Langevin equation can be derived in a generic way, 
so-called the Mori theory [12, [12] , from the microscopic 
equation of motion. The starting microscopic equations 
of motion are Liouville or Hisenberg equation for a clas- 
sical or a quantum system, respectively. They read 

^^A,(t)^{Mt),H}pB, (5) 
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and 



d_ 

di 



A,{t) = (l/ih)[A,{t),H], 



(6) 



respectively. Here, Aj{t) are arbitrary slowly varying 
variables (slow variables), H a microscopic Hamiltonian, 
and, {•, 'Ipb and [•, •] represent Poisson brackets and com- 
mutation relations, respectively. Equation ([5]) (or Eq. 
([6])) can be divided into two parts: one is composed of 
only the slow variables and describes their slow motion, 
while the other involves fast motions due to the micro- 
scopic processes. Thus, we have [IS, HI] 



5Ak 



+ e,{t), (7) 



with j3 being the inverse temperature. Here, Vj{A(t)), 
which is called a streaming term, gives the time-reversible 
process, while Ljk and H{A{t)) are bare kinetic coef- 
ficients and a thermodynamic potential for Aj^ respec- 
tively. The first and second terms in Eq. (O are the slow 
motions and nonlinear in Aj , whereas the last term is the 
fast motions and treated as a stochastic variable obeying 
the fluctuation-dissipation relation 



{0j{t)9kit );a> = 2Ljk{a)S{t-t 



(8) 



Here, 



; a) represents the conditional average in which 



Aj is fixed at aj . We stress that this relation is not given 
by hand, but naturally obtained from the decomposition 
process |22| . 

Equation ([7]) is called the generalized Langevin equa- 
tion, which has been widely used in the phase transition 
dynamics [13, [HI- Even for the QCD CP, we may use the 
generalized Langevin equation, because only the time- 
scale separation is assumed in the Mori theory. Further- 
more, we note that, by the time-scale separation, trans- 
port coefficients arises . 



C. Dynamic RG 

The general dynamic RG transformation usually con- 
sists of two procedures, i.e., coarse graining and rescaling 
as in the static RG transformation [21,,^]. However, as 
is shown in [13, [13, [l^j , we can omit the rescaling, if we 
are interested in only the critical exponents of transport 
coefficients, although the relevant fixed point seems to be 
absent in such a simplified RG transformation [33 |. 

The Langevin equation is an infrared effective theory 
and inherently has a ultraviolet cutoff Aq, which should 
satisfy the following inequality 



C ^ < Ao < a" 



(9) 



Here, ^ and a are the correlation length and the char- 
acteristic microscopic length, respectively. Then, the 
Langevin equation is coarse grained by averaging over 



the large-wavenumber components of the slow variables 
Aj{t) in the infinitesimal wavenumber shell. 



A - (5A < A; < A, 



(10) 



for Eq. ([7]). Here, A starts from the initial value Aq and 
is lowered up to A <C In this way, we infinitesimally 
make coarse graining of the Langevin equation. Because 
the coarse-graining procedure is infinitesimal, we need 
not the rescaling . Inspecting the form of the coarse- 
grained Langevin equation, we can obtain the RG equa- 
tion for the transport coefficients. 



D. Contrast with the static RG 

Here, we first stress that the concept of the dynamic 
universality class is not so universal contrary to its name. 
Then, the class of the QCD CP may not be the same as 
of the liquid-gas CP or the model H although it is 
conjectured by [H, [l^]- To see this, let us contrast the 
difference between the static RG with the dynamic one. 

An important point is that the respective infrared ef- 
fective theories are different; in the static case, the in- 
frared effective theory is the thermodynamic potential 
(or so-called Landau free energy), the nature of which 
turns out to be governed by only the space dimension 
and the symmetry among the order parameters but not 
by the details of the dynamics, and hence the concept 
of the universality class makes sense for the static case. 
In contrast, for the dynamic RG, the infrared effective 
theory is the nonlinear Langevin equation. 

Here, the important difference arises; the relevant vari- 
ables for the Langevin equation is not only the order pa- 
rameters but also conserved densities, and its nonlinear 
couplings can not be determined by only the symmetry 
in general. Consequently, the dynamic universality class 
is not so universal compared to the static one. Specif- 
ically, the nonlinear couplings, namely, the streaming 
terms Vj{A), are generally given by the Poisson brack- 
ets (commutation relations) am ong the slow variables in 
the classical (quantum) system [2l|; 



v,iA) = ^[0,.,(A)|£ - /3-i_i_g^.,(A)], 

A: 



(11) 



where 



Qjk{A) = {{A,,Ak}pB;A) or {[A,,Ak]/{in);A). (12) 

The important point is that the Poisson-bracket relations 
depend on the microscopic expressions of the variables. 



^ The model H |25l . l27l | is the minimal-dynamic model for a CP that 
its relevant modes are given as the nonrelativistic-hydrodynamic 
modes. The liquid-gas CP belongs to the dynamic universality 
class of the model H 
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This fact leads to an important consequence that the 
dynamic universahty class of the QCD CP may not be the 
same as of the liquid-gas CP or the model H. Actually, in 
the model H, the Poisson-bracket relations are calculated 
with the non-relativistic relations |2ll u3i ■ 



III. THE NONLINEAR LANGEVIN EQUATION 
FOR THE QCD CP 

A. Slow variables 

We first identify the slow variables near the QCD CP, 
which consist of soft modes and conserved densities. For 
the QCD CP, the soft modes are nothing but the long- 
wavelength fluctuations of the conserved densities, i.e., 
the baryon number density n and the energy and momen- 
tum T^i, [ij]. Thus, we see that the slow variables 
near the QCD CP are given by only the fluctuations of 
the conserved quantities: 



Aj = {6n,6e = {ST°"),Sr = {ST°')}. 



(13) 



Because the slow dynamics of the conserved quantities 
is basically given by hydrodynamics, we find that the 
system near the QCD CP is described as a relativistic 
critical fiuid. In other words, the relevant modes are 
given as the relativistic-hydrodynamic modes. This is 
the basic observation for our construction of the nonlin- 
ear Langevin equation for the QCD CP. More specifi- 
cally, the hydrodynamic modes are the thermal and vis- 
cous diffusion modes, and the sound propagating mode. 
The thermal mode is the entropy fiuctuation inducing 
the density and energy fluctuations, whereas the viscous 
and the sound modes are the transverse and longitudinal 
components of the momentum fluctuations, respectively. 

Now, we note that not all fluctuations are enhanced 
near the CP. Therefore, we can neglect nonlinearity of 
fluctuations that is not enhanced, if such fluctuations are 
identified. Then, let us identify the non-enhanced fluctu- 
ations by the hydrodynamics. The usual hydrodynamics 
with static scaling laws is useful to see the such tendency 
of the slow variables. Since the result turns out to be 
independent of the choice of the frame, which defines the 
local rest frame [1^, let us take the energy frame, namely 
the Landau equation [23] , which is given by the following 
conservation laws: 



0, 

= 0, 



(14) 
(15) 



where N'^ and T^^^ are the particle current and the 
energy-momentum tensor, respectively. Those are given 

as 



(16) 
(17) 



where h = e -I- P is the enthalpy density with e and 
P being the energy density and the pressure. Also, 



u^^ = (7, 71') are the fluid four velocity, with 7 being the 
Lorentz factor, and the dissipative terms, ly'^ and t'^" , 
are given by 



+CoA^-{d^-u), (19) 



where Aq, r]Q and Co are the bare thermal conductivity, 
the bare share and bulk viscosities, respectively. A'^^ = 
gf^f yMyi^ ig the projection onto the space-like vector 
and 5^ = A^'^di, is the space-like derivative. 

In the hydrodynamic regime, fc^ ^ 1, the hydro- 
dynamic modes is analyzed by the linearized equation, 
which is given by 



dSn 

'dt 

dSe 

'dt 
dSJ 



— — ——ncV ■ Sv + 



nrTr 



T 



V^<5^ , (20) 



(21) 



=-V(<5P) + (Co + 3r?o)V(V • 6v) (22) 



where the symbols with a preflx 6 denote the fluctua- 
tions from their equilibrium values, which are denoted 
by a suffix c ^. Hereafter, variables with the suffix and 
the prefix respectively denote the equilibrium values and 
fluctuations. As relativistic effects, we see that dissipa- 
tive effects appear in Eq. (|20p while vanish in Eq. (|2ip . 
because we have chosen the energy frame. We note that 
the relativistic effect for the particle frame appears in a 
different form '23']. 

By Eqs. ([20l) - ((23|) and static scaling laws, the tendency 
of the hydrodynamic modes is analyzed in the critical 
regime, fc^ ^ 1. We have studied the such behavior in the 
previous paper [23j and shown that the thermal mode is 
enhanced , whereas the sound mode is suppressed and the 
viscous mode is not enhanced nor suppressed. Recalling 
the relation between the hydrodynamic modes and the 
slow variables, we have the result that Sn and Se are 
enhanced, while 6 J is not near the QCD CP. 

We note that nonlinear couplings among these modes, 
which is not included in usual hydrodynamics, become 
significant in the critical regime. We will take them into 
account in the nonlinear Langevin equation, except for 
the fluctuation of the momenta SJ, the nonlinear term 
of which will be neglected even in the critical regime. 



* Here, we have slightly rewritten the form of the linearized equa- 
tions, 1 I2OII - II23I I. from those in [23| by the thermodynamic rela- 
tions Se = T5{ns) + /iSn and SP = nsST + nSfi, where s is the 
entropy per particle. 
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B. Thermodynamic potential for the slow variables 

Next, we construct the thermodynamic potential 
H{6n, Se, SJ) for the slow variables. 

Since the momentum density fluctuation is not en- 
hanced near the QCD CP, we can neglect its coupling 
with Sn and Se, and may adopt a Gaussian form for the 
momentum density part of the potential. Thus, we have 
H{Sn, 5e, J) ^ Hne{Sn, 5e) + Hj{5J), with 



Hj{SJ) 



1 

2h, 



■SJ' 



(23) 



In contrast to SJ, Sn and Se are significantly en- 
hanced near the QCD CP, the thermodynamic potential 
Hne{Sn, Se) should contain higher order terms of them. 
In fact, Hne{Sn,Se) is the quantity that determines the 
static property of the system and the QCD CP belongs 
to the same static universality class as the 3d Ising class, 
namely, Z2. Therefore, we may construct H„e{Sn,Se) 
with the thermodynamic potential for the 3d Ising sys- 
tem [33, which reads 



l3Hising{ip,m) 



-1-70-0^ »7i + T^m"^ 



hip- 



(24) 



Here, and m are the spin density and the exchange 
energy density, respectively, tq, Kq, uq, 70 and Cq denote 
the static parameters, while h and r the applied magnetic 
field and the reduced temperature, respectively. Then, 
we assume that the thermodynamic potential takes the 
following form 



H{Sn,Se,SJ) = i?ising(V') ™) 



1 

2h, 



-SJ' 



(25) 



provided that the mapping between (tp, m) and (Sn, Se) 
is given. 

The general mapping relation between a grand canon- 
ical ensemble in Z2 and the 3d Ising system is known in 
condensed matter physics (2Qj . which are summarized as 
follows. First, we assume the following linear relation be- 
tween the deviations of the respective intensive variables 
from those at the critical points, which should be valid 
near the CP. ^: 



Sh —aiS{fi/T) + a2ST/Tc 
St ^l3iS{fi/T) + p2ST/T„ 



(26) 
(27) 



where ai, a2, Pi and (32 are constants and assumed to 
be regular at the CP. We note that ai, a2, Pi and j32 
need not to be determined for the critical divergence of 



transport coefficients, since they have no singularities at 
the CP. Although one could use Eqs. (US]) and (gT]) for the 
mapping, it turns out to be inconvenient for the analysis 
by a Langevin equation. To translate these equations to 
more convenient ones, we assume the following relation 



iPSh + mSr = T'^STSe + S{n/T)Sn, 



(28) 



which is actually derived by considering a change of the 
microscopic distributions due to small deviations of the 
external parameters in both systems. From the relations 
Eqs. (I26l) - (|28l) . we arrive at the convenient mapping re- 
lation as follows, 



Sn —aiip + Pim, 
T^^Se —a2t]] + /32W. 



(29) 
(30) 

With this mapping, Eq. (|25p now gives the thermody- 
namic potential for the QCD CP. We remark that we only 
map the static quantities although the dynamic ones are 
studied. For later uses, we introduce fluctuations of the 
intensive variables as 



Se - 
Si(3H) 
TJ Sn ' 



(31) 
(32) 



This relation comes from the fact that, in the grand 
canonical distribution Pg^a cx: exp[(l/r)e + (ii/T)n], e 



and n are respectively conjugate to 1/T and /i/T 
We also introduce the fluid velocity fluctuation as in the 
non-relativistic case: 



Sv 



S_H_ 
SJ' 



(33) 



We note that the static parameters in Eq. (|24|) has the 
ultraviolet cutoff dependence in the region < A. Let 
us write the static parameters as r(A), K{A), u{A), 7(A) 
and C(A) to make their A dependence explicit. These 
variables have the following asymptotic behaviors [13, [H, 
M: 



riA)r 

K{K)r 
u{K)r 

7(A)^ 
C(A). 



A2-", 

A-", 

A'-2^, 

^(e+a/i.)/2-,, 

A-"/", 



(34) 
(35) 
(36) 
(37) 
(38) 



where e = 4 — d with d being the space dimension, while 
a, V and 77 are the usual static critical exponents. Noting 
that 77 is of order and very small, we neglect r\ and set 
— 1, hereafter. 



C. Streaming terms and bare kinetic coefficients 



® Recall that the static scaling laws are expressed by the deviations 
of the intensive variables from those at the CP. 



In this subsection, we determine the forms of the 
streaming terms, d„, v^^. and vj. We can nicely deter- 
mine the flrst two terms from the continuity equations 



6 



because Sn and Se are the conserved densities. From the 
continuity equations, we can write w„ and as diver- 
gences of reversible currents, which read 



(39) 
(40) 



with j„ and being the reversible currents of the 
number and energy density, respectively. Here, 7 is 
the Lorentz factor of the fluid-velocity fluctuation, n = 
Tie + Sn and e = ec + 5e. As the reference frame, we have 
chosen the rest frame of the equilibrium state, and then 
the back ground fluid velocity vanishes. Furthermore, 
We may set 7 ~ 1, because the fluid- velocity fluctuation 
is given by Sv = h~^SJ that is not enhanced. Therefore, 
we write the streaming terms, Vn and v^, as 



Vn= -V • (nSv), 

Ve= -V ■{{e + Pc)5v), 



(41) 
(42) 



where we neglect the pressure fluctuation because it is 
not enhanced near the CP [23j . 

Now, we note that the determination of vj is not sim- 
ple. Although the continuity equation tells us that vj is 
the divergence of the reversible-stress tensor, the determi- 
nation of the reversible-stress tensor is not trivial. How- 
ever, we can determine it from the potential condition, 
which is a general condition for the streaming terms [20| . 
The potential (or divergence) condition 20, 21] reads 



dr 



j—7i,e,J 



6A, 



dr ^ 



(43) 



We remark that this condition can be derived from the 
definition of streaming terms : 



v,iAit))^{A,;A{t)), 



(44) 



where Aj = iCAj (0) is the microscopic time derivative of 
Aj and iC is the Liouville operator. In a continuum sys- 
tem, the right-hand side of Eq. (|43l) vanishes in general 
[2(]j . Thus, the potential condition is reduced to 



0, 



(45) 



j—n.e^J 



where Vj is only the unknown quantity because we have 
already determined Vn, Ve and H{Sn, Se, SJ). Using Eqs. 
([M]). (HT|) . dUl) and (gS]), we obtain 

/SH SH 
dr[n\7— + (e + Pc)V— + vj] ■ (3Sv = 0. (46) 
on de 

Since this condition should be satisfied for an arbitrary 
fluid-velocity fluctuation, we have 

vj = -nV^ - (e + Pc)^^. (47) 
on oe 

Next, let us determine the kinetic coefficients from 
the relativistic hydrodynamic equation, Eqs. (ITH) - (ITO1) . 



From Eqs. (US]), ^9^, §2} and dM]), we can read the 
kinetic coefficients Ljk for small Sv as 



(48) 



-TMjd.dj + (Co + (1 - 2/d)77o)a,9,], (49) 

and that the other coefficients are zero. We note that 
Lee is absent due to the choice of the energy frame. 

Now, we have determined all the necessary terms, and 
then can write down the nonlinear Langevin equation for 
the QCD CP as 



dSn „ , , , , Si^H) 

^ =-V-iie + Pe)Sv), 

dSJ SH , , SH 



(50) 
(51) 

(52) 



where 0„ and 9j are the noise terms and satisfy the 
fluctuation-dissipation relations 



,ir,t)0n{r ,t )) ^-2Ao 



TlrTr 



V'^S{r ~ r 



(53) 



xS{t - 1 ), 

+{Co + (1 - 2/d)f]o}d'd^] 
xS{r -r')S{t^t). (54) 

Let us write the transport coefficients as A (A), 77(A) 
and C(A) to make their cutoff dependence in the region, 
< A. The critical behaviors of the transport coef- 
ficients are determined from their asymptotic behaviors 
near the relevant fixed point as A is lowered. 

Here, we compare the Langevin equation Eqs. ([50)1 - 
(1521) with that for the liquid-gas CP ^ 



9Sn „ , ^ , 
— =-V . inSv), 



(55) 



^ =-V • ((e + Pe)Sv) + XoTeV'—- + Oe, (56) 
ot Se 



dSJp 
dt 



-nV 
-Ljj- 



JH^s 



Sn 
S{(3H^' 

SJo 



(e + Pc)V 



~sr 



+ 0j, 



(57) 



where SJp = PcSv, p and H^^ are the non- relativistic 
momentum density, the mass density and the thermody- 
namic potential for liquid-gas CP, respectively: 

H'^iSn, Se, SJ) = iJisi„g(^, m) + -^SJ^. (58) 
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We see that the streaming terms have the same forms but 
the dissipative ones are totally different between the rela- 
tivistic and non-relativistic cases. The difference also ap- 
pears the relation between the momentum and the fluid- 
velocity fluctuation. Therefore, one may naturally expect 
some novel characteristics in the relativistic case that is 
absent in the non-relativistic case [2a|. 



IV. THE TRANSPORT COEFFICIENTS BY 
DYNAMIC RG 

We here present an analysis of transport coefficients 
by the dynamic RG. A detailed derivation of the RG 
equations is given in the Appendixes. 

First, we rewrite Eqs. ([50l) - ([52]) as the equation for 
-0 and m to conform the hydrodynamic variables, 5n and 
(5e, to the Ising variables, "0 and m. Noting that we can 
set a2 = in the mapping relations, Eqs. ([^ and ([SHI), 
without loss of generality [1^ , we have 



d'4' 



6m) 

' <50 



-h-^V ■ (tpSJ) 



— =-^2"iV • 6 J - h-^V ■ [mSJ], 



dm 
'dt 
d{SJ) 
dt 



(59) 
(60) 



dip dm 
oip dm 



(61) 



with C^, = a^'^{nch'^^ — I3i(32^) and Cj = a^^{nc 
(3ihc). Here, we note that we could rewrite the potential, 
Eq. ([24|) . for and m as that for Sn and Se to conform 
the variables; the choice is a matter of preference. 

In the dynamic RG transformation, we average over 
the short wavelength components in the shell, A — 6A < 
k < A, for the Langevin equation. For this task, we 
must perturbatively solve the equation about them, by 
rewriting it as a self-consistent equation [2l|. Although 
an explicit derivation of the self-consistent equation for 
the QCD CP is first made in this paper, we leave the de- 
tails of the derivation to Appendix|3 because the general 
procedure of the derivation is standard and given in the 
textbook (21(]. Here, we present only a few basic equa- 
tions of the dynamic RG for the QCD CP. Now, as is 
shown in Appendix lAl Eqs. (|59 l) - (j6l]) can be reduced to 
the following form; 



■0(fc,a;) 
m{k, uj) 
ySJ\\{k,uj) 



I +G"{k,uj)V{k,uj), (62) 
ySJ°ik,u;)^ 



and 



where J|| (A;) = k- j{k) and J±{k) = J{k) — J|| (fc) are the 
longitudinal and transverse components of the momen- 
tum. Here, and G*[ are the bare propagators, which 
are given by Eqs. (IA50p and (|A5ip - (|A53I) . whereas 
V and V the nonlinear couplings, coming from the 
streaming terms and given by Eqs. (fMO)) - (|X34l) and 

(IA59p . Also, 0°, 5J^i^ and 5J^ are the bare variables, 
which are the solutions without the nonlinear terms. It- 
erating the self-consistent equations (|62|) and (|63|) . we 
can obtain a perturbative expansion of the nonlinear cou- 
plings and have a coarse-grained Langevin equation. 

Now, we note that the variables, V',Jj_ and J||, are re- 
spectively correspond to the thermal, viscous and sound 
modes ^ (see the propagators (|A5ip - (|A53|) .'). Therefore, 
the first and third rows of Eq. (15^ respectively denote 
the equations of motion for the thermal and sound modes, 
while Eq. ((63|) for the viscous mode. We stress that the 
sound mode is neglected in the model H, although it is 
essential for the renormalization of the bulk viscosity. 

Here, we make the coarse graining to the second order 
in the nonlinear couplings, V and Vj_^^ (see Fig. [3] for 
an example.). Inspecting the coarse-grained equation for 
tp (see Eq. ()B5|) for the detail) and, we have the RG 
equation for the thermal conductivity: 



A 



dX{A) 
OA 



r/(A)A(A), 



(64) 



/(A) = TcK4/{D^r]{A)X{A)A^), is the surface area of 
a unit sphere in 4 dimensions divided by (Stt)**, = 
(ricTc/aihc)^. Here, we have introduced /(A) for conve- 
nience sake. Similarly, from the coarse-grained equations 
for Jx and J|| , we obtain the RG equations for the shear 
and bulk viscosities 



-A 



dC{A) 
dA 



--Aj^{A)X-\A)A' 



■e-4 



(65) 
(66) 



where 7(A) is a static parameter in the thermody- 
namic potential (see Eqs. ([24| and ([37| ). and A = 
h'^^Ki/ {j32D^). Furthermore, differentiating /(A) about 
A, we also have the RG equation for it: 

-A^-/(A)(e-g/(A)). (67) 

Now, we note that Eqs. (IMl) . and ([57)) are identi- 
cal to those for th e liq uid-gas CP except for unimportant 
constants in ffA)|20l [26|. Equation is also equiv- 
alent to the RG equation of the bulk viscosity for the 
liquid-gas CP in the limit cj ^ [13, HI]. Therefore, ar- 
guments about the RG equations and results from those 



5J^{k,oj) = 5J^(k,uj) + G\V^^^{k,uj), (63) 



^ Although rh would be a linear combination of the thermal and 
sound modes, we need not to consider rh for a following analysis. 
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are the same as for the hquid-gas CP. Then, we provide 
only essential arguments and results in the following part, 

and leave the detail to [13, HE im, HI. 

Now, we identify the relevant-fixed point as the fol- 
lowing [2^. Because, at a fixed point, parameters 
are invariant about the RG transformation, we set the 
left-hand side of Eq. (1571) as 0. Then, as a fixed-point 
value of /(A) which is denoted by /*, we have /* = 
and /* = (24/19)e. Therefore, we have the two fixed 
points and the relevant one is specified by /* = (24/19)e. 
Although the relevant point seems to be absent in Eqs. 
(|64| . ([65]) and ([66|). the reason is due to the simplified 
RG transformation as mentioned in the earlier section, 
and this is just a apparent problem [sills^. 

Substituting /* = (24/19)e into Eqs. ([Ml), §B and 
(|66p , we have the asymptotic behaviors near the relevant- 
fixed point: 



A(A) -A^iS 

77(A) --A^T9 

C(A) ^A-(4 



7) 



(68) 
(69) 
(70) 



Here, in the derivation of Eq. ([70]). we have used the 
asymptotic behavior of 7(A), Eq. ([57| . Decreasing the 
cutoff to the region A ^ £~^, we can replace A with i^"^ 
in the asymptotic behaviors |20l . Isij : 



Ar 



In three dimensions, we find 



tO.95 



(71) 
(72) 
(73) 



(74) 
(75) 
(76) 



We can also read the dynamic critical exponents from 
Eqs. ((7T|) - ([75)) . A dynamic critical exponent, denoted 
by z, generically parametrizes the decay rate r(A;) at the 
wavenumber k = as r(£~^) ~ shown in 

Appendix the decay rates for the three modes at k 
are given by 



rthGrmal(fc) ^^^^^(rR + k^)D^, 
rviscous(fc) =?7Rfc^/?-^\ 
r.ound(fc) =(Cr + 2(1 - l/d)7^n)k''h-\ 

Thus, wc find the dynamic critical exponents as 



18 

^^thcrmal — ^ 19^' 



-^VISCOUS 



-S^sound 



19 



-(2 



18 
19^ 



(77) 
(78) 
(79) 



(80) 
(81) 
(82) 



In three dimensions, the dynamic critical exponents are 
given by 



^thermal '^3, 
■^^viscous 
-^^sound 



-0.8. 



(83) 
(84) 
(85) 



We see that the thermal and viscous modes exhibit criti- 
cal slowing down, while the sound mode critical speeding 
up. 

Why do not the relativistic effects appear in the RG 
equations? The reason is that the nonlinear terms in 
the dissipative terms generally renormalize only static 
parameters, up to order [2l|, |3l|. Furthermore, the 
difference in the relation between the momentum and 
the fluid velocity is only unimportant constants, i.e., the 
enthalpy density h and the mass density p. Then, the 
RG equations are essentially the same as for the non- 
relativistic case. 



V. SUMMARY AND CONCLUDING REMARKS 

We have studied the critical behaviors of the trans- 
port coefficients and the dynamic critical exponents at 
the QCD critical point (CP) by dynamic renormaliza- 
tion group (RG). For this purpose, we have constructed 
the nonlinear Langevin equation near the QCD CP for 
the first time. Our construction is based on the gener- 
alized Langevin theory, by Mori [l^, [2^ , and the rela- 
tivistic hydrodynamics; instead of a naive construction 
method [2l| , we have determined the streaming terms by 
the relativistic hydrodynamics and the potential condi- 
tion that gives a constraint to these terms. The resulting 
equation is given by Eqs. ([50| - ([52|) . Although there are 
some attempts to make a one-to-one mapping between 
QCD CP and Ising CP [il,!!^, we have shown that it is 
not necessary to specify such the mapping for the critical 
exponents, as for the liquid-gas CP 20]. 

We have shown that the bulk viscosity and the thermal 
conductivity strongly diverge at the QCD CP. Also, we 
have found that the thermal and viscous diffusion modes 
exhibit critical slowing down with the dynamic critical 
exponents zthcrmai ^ 3 and ^viscous ^ 2, respectively. In 
contrast, the sound propagating mode critical speeding 
up with the negative exponent Zsound ^ —0.8. 

We now compare our result about the bulk viscosity to 
that in 16]. Although a divergent behavior of the bulk 
viscosity is the sarne, the critical exponents is different in 
the two cases. In [1^1, the critical exponent is estimated 
to be about 0.2 and the divergence is weak contrary to 
our result. We also note that the study by the relativistic 
Boltzmann equation 19] gives only the bare bulk viscos- 
ity. 

We note that the bulk viscosity and the thermal con- 
ductivity are usually neglected in heavy ion physics, how- 
ever they become much more important than the shear 
viscosity near the QCD CP. Furthermore, the description 
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for the created matter as a perfect fluid is not valid near 
the QCD CP at all due to the strong divergence of the 
bulk viscosity. 

As the argument about the dynamic universality class 
[3 13 : we have shown, from an explicit calculation, that 
the QCD CP has the same critical behaviors as the liquid- 
gas CP has. The argument assumes the insignificance of 
the relativity for the critical dynamics by the slowness of 
the diffusion processes. However, we have shown that the 
genuine reason for the insignificance originates from the 
small fluctuation of the momentum density; the critical 
dynamics is essentially governed by the streaming terms, 
which are modified by the relativistic effect through only 
a Lorentz factor of the fluid velocity fluctuation. How- 
ever, the fluid velocity fluctuation, which is proportional 
to the momentum, is not enhanced near the CP. Thus, 
the relativistic effects do not affect the critical dynamics 
near the QCD CP. We stress that the sound mode exhibit 
critical speeding up, and then the sound diffusion is fast 
near the QCD CP. Therefore, the basis of the conjecture 
would be true for the thermal and viscous modes, but not 
for the sound mode. We also note that the model H 27], 
which is the mmimaZ-dynamic model for the dynamics 
near the liquid-gas CP, can not describe the critical be- 
havior of the bulk viscosity because it does not contain 
the sound mode. 

Wc note that our Langevin equation must satisfy usual 
fluctuation-dissipation relations, Eqs. ([55)) and for 
the consistency with the linearized Landau equation ^, 
although a relativistic Brownian motion seems not to 
satisfy the usual relations [sj]. Moreover, our Langevin 
equation seems to violate the causality, since the dissi- 
pative terms are determined from the Landau equation. 
However, the Israel-Stewart equation [ssj . in which the 
causality problem is formally resolved, gives the same re- 
sult as the Landau equation gives in long-wavelength re- 
gion [IJ . Therefore, our determination from the Landau 
equation must suffice. Furthermore, we note that short- 
wavelength components in the region, k > where a 
is a characteristic microscopic length, would violate the 
causality. Therefore, we can exclude such the illegal com- 
ponents from the theory by the cutoff, Aq <C a^^. We 
stress that all infrared effective theories inevitably have 
the ultraviolet cutoff; naturally, relativistic hydrodynam- 
ics also has it. 

Also, we note a frame dependence of the results. As 
a hydrodynamic equation, we used only the equation in 
the energy frame. Does the results change if an equa- 
tion in the particle frame is used? Although the frame 
dependence can appear in only dissipative terms, the crit- 
ical dynamics is essentially determined by the streaming 
terms. Therefore, the results would not change for the 
particle frame, if an equation in the frame is correct. 



However, in practice, the Eckart equation has a patho- 
logical behavior [3^. Namely, fluctuations do not relax, 
and therefore we cannot use the Eckart equation. 

Recently, some authors have suggested the existence of 
other critical points in higher density region of the QCD 
phase diagram where the color superconductivity is taken 
into account 0, Q • It would be interesting to study the 
critical dynamics near such a new QCD CP using the 
dynamic RG theory, as an extension of the present work. 
For this purpose, however, we must firstly specify the soft 
modes and construct the nonlinear Langevin equation. If 
the soft modes are different from the conserved densities, 
which is the case when the the diquark fluctuations are 
relevant [l, [s^ , the construction based on the relativistic 
hydrodynamics done in the present work does not work, 
and we must directly recourse to Eq. (jlip to identify the 
streaming terms. 
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Appendix A: Rewriting the nonlinear Langevin 
equation as a self-consistent equation 



Here, we rewrite the Langevin equation, Eqs. (|59t - (|6T|) 
as a self-consistent equation. First, we make a Fourier 
transformation as the following 



'4>{k,uj) 



dtd'^re''^'~'^-^ij{r,t). (Al) 



Then, we have 



If our nonlinear Langevin equation is linearized, the linearized 
equation must give the same result as the Landau equation gives. 



— ILJipik^Lj) — — L^ik ■ oJ — a. Lnn — ; 

-h-\k- f {'4^{q)SJ{k-q))+a^'0llA2) 

—iwrh{k,u}) =—(32^ik ■ 6 J 

-h-^ik- [ {m{q)SJ{k-q)), (A3) 
Jqn 

—iuj5J{k,uj) ——Cjik^ — j32^hcik- — 
oip om 

(A4) 



-{T,h,)~^Ljj-6J + e, 
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Note that the quantities with tilde in Eq. (IA2|) - (IA4|) are 
Fourier transformed, like Eq. (|Aip . and we have abbre- 
viated the nonlinear terms such as 



^{q)5J{k-q) = / ——^4,{q^n)5J{k-q,uj-n). 

(A5) 

We now decompose Eq. (IA4p into the longitudinal and 
the transverse components: 



SH 



SH 



iuiSJu ——iCjk— i/So ^hck- — 

" dip dm 



-I (k- q) 

Jqil 

X [it(«)^(^ -q) + Y^('7)"i(fc - q)] 
dip dm 

-{T,H,y^k-Lj.j{k)-5J +~e\\, (A6) 
~iuj5J^= —i / V^{k) ■ q 

dtp dm 
-{T,h,)-^V^{k) ■ Ljj{k) ■ 5J 
+9±, (A7) 



where we have introduced a projection operator as 



and 



irA_{k)),j^d,j~hk,/k^ 



SJ\i{k)^k-SJ{k), 
6J±{k)=^V±{k)-SJ{k), 
0\\{k)=k-e{k), 

e^{k)^v±{k)-e{k). 



(A8) 



(A9) 
(AlO) 
(All) 
(A12) 



Because the streaming terms in Eqs. (|A6p and (IA7[) are 
too complicated for our purpose, let us retain only the 
terms that yield dominant contributions for the transport 
coefHcients. We note that only such terms suffice for 
obtaining the critical exponents. From the relations [20j 



dV(m(r)m(0)}~ 



(A13) 
(A14) 



we expect tp yields stronger singularity than m. There- 
fore, we only retain the term that are of the second order 



in ip. Namely, we reduce the streaming terms as 

iCjk^ + i/S^^hck^ 
dip dm 

f S H S H 

+1 (fc • q)[—{q)i!{k -q) + ^{q)m{k - q)] 
' dip dm 



qU 



Tc[iCjkxo\k)iP + ip-^hckC^^ih 
+iP2^hck^o ip{q)ii}(k - q)], 

J qO. 



(A15) 



f 6H SH 

« / ^i(fe) • Q[^i<l)i'ik -q) + ^{q)rh{k ~ q)\ 
Jqn dip dm 

^iT,r^{k)- ( qx^^qMqmk^q), (A16) 
Jqn 

where Xo ^(fc) — fQ + k"^. Notice that we have set Kq = 1, 
as mentioned in the text. 

Next, we consider the dissipative terms. The impor- 
tant point is that the nonlinear terms in dissipative terms 
renormalize only static parameters in a thermodynamic 
potential to second order in e, generally [2l|,|3l|. There- 
fore, we can take into account nonlinear terms in the 
dissipative terms with the results of static RG, Eq. ([M)) - 
(l38l) . and effectively neglect it in the Langevin equation. 
Then, we reduce the LnnS{/3H)/Sip as 

L„„(fe)^(fe,c.)^Aofc2xoi(^^^ (A17) 

In contrast, the dissipative terms of dJ are originally 
linear and then directly read 

k ■ Ljj{k) ■ 5J{k,u)= TjCo + 2(1 - l/d)7?o] 

xfc25J||(fc,w), (A18) 
V^{k) ■ Ljj{k) ■ dJ{k,uj)= T^TjoP6J^{k,uj). (A19) 

Collecting the above results, we arrive at the reduced 
nonlinear Langevin equation as given by 

— iujip ——ikC^dJ\\ 



^h-^ik- / iP(q)dJ{k~q) 
Jqn 

^\ok'^D^XQ'^{k)ip + a^^On, 



—iujm ——B^ ^ikdJw 
Jqn 

-JW(5J|| ~Tc[—ikxQ^ik)Cjip — ikC^^ I32hcm 

-ikl3^^hc^o / '4'iq)4>ik - q)] 
Jqn 

-k'^v\^h~'^dj\\ +e\\, 



(A20) 
(A21) 



(A22) 



-iLodJ^ =-iT,v^{k) ■ / gxo ('z)V'(g)V'(fc - q) 

Jqn 

-k^f]oh-'^dJ^+e^_, (A23) 
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where 



D, 



iy'o= [Co + 2(l-l/d)r;o]. 



(A24) 
(A25) 



This is the basic equation for the dynamics near the QCD 
CP, which is first written down, and a main result of this 
paper. 



We can compactly rewrite the basic equation in a ma- 
trix form: 



M{k,uj) Th{k,uj) 
VW||(fc,^)y 



where 



V{k,Lu) + e{k,Lj), (A26) 



M{k,uj) = \ 

ikxo^ik)CjTc ikCg^ f32hcTc-iuj + k^iy^^h-^ 





-iuj 



tkCip 
ik(i^^ 



(A27) 



6»(fe,w) 



' a-[^e{k,uj)^ 




V^^±{k,uj) + V^^\\{k,uj) 

V{k,ijj)= ( Vmm^{k.,Uj) + Vmni\\{k.,Uj) 

V||^^(A;,w) 



(A28) 



(A29) 



and 



V^^\\{k,uj 

Vmm\\{k,i^ 
V||^^(fe, W 



= -h-^k ■ / i^{q)SJ±{k - q), (A30) 



XiP{q)SJi\{k~q), 



f k-{k-q)/\k-q\ 



-h,. ^ik 



(A31) 

m{q)dJ^{k-q),{A32) 



qfl 



-h-^ik / m{q)5J\\{k-q), (A33) 
-ikT,l3-^K-fo U{q)i,{k-q). (A34) 

J qil 



Since Eq. (|A23p is decoupled from the other equations 
at linear level, we do not rewrite it as the matrix form. 

Next, we calculate the bare propagator G'^{k,uj) — 
A4~^{k,u!). The inverse matrix is given as the trans- 
posed cofactor matrix divided by detTW. The determi- 
nant reads 

detM ={~i^f + {-ii^fk^{\D^Xa^{k) + i^oK^) 
-icj^iCo^h^Tc + XaHk)C^CjTc) 
+k^XoXo^D^Co^h,Tc 

-i^k^\oX^^D4k)vih-\ (A35) 

Here, in the coefficient of —iajk"^, taking into account the 
behaviors after renormalization [26l l33j | which are given 



r 



as 



-0.2 



XB'ik)-r^ + k', 



(A36) 
(A37) 



we neglect Xo ^{k)C^CjTc by comparing with Cq ^hf-T^ 
. Then, we can factorize the determinant as 

detA^- {~iuj + Xo{k)xoHk)) 

X {—ioj + ikcs + ^i^o^c^k^) 
1 



X {-iuj - ikCs + 2 o^<^ 



(A38) 

in the long-wavelength region. Here, we have defined 



(A39) 

c;=Ua'ncic- (A40) 
The diagonal components of the cofactor matrix m reads 

1 

— I 
2 



Ao(fc)= Xok-^D^, 

Cq hcTc- 



mil ^{—iuj + ikcs + —i^gh^ ^k"^) 



(A41) 



=(-icj)^ - iujk'^ XQXo{k)D^iyQh^ ^ 

+k^X.o\k)C^CjT, 

+k'^Xaxo{k)D^iyQh^^, 
={-iuj){-iuj + XQ{k)xQ^{k)), 

and the ofF-diagonal components are given by 

mi2=ex-Hk)P2'CjT,, 
mi3= -kojXQ^{k)CjTc, 
m2i= -k'^CQ^hcC'^l32Tc, 
77123= -ikC^^KP^Tci-iuj + kHoXo{k)D^), (A47) 
m3i= k^Ca^h,C,pP2Tc, (A48) 
77732= -ikp2\-^i^ + Xo{k)xo\k)). (A49) 



(A42) 
(A43) 

(A44) 
(A45) 
(A46) 
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Here, we neglect the ofF-diagonal components because 
they would not yield dominant contributions to the trans- 
port coefficients. Then, we obtain the bare propagator 
as 



Using the fluctuation dissipation relation Eq. ((53)) . we find 



GO(fc,c.) 
G°(fc,w) = I G^(fc,a;) 

G°ik,uj)^ 



with 



G°(fc,a;) 



G||(fc,a.)^- 



1 



-iui + ikcs + ^k"^ 
1 



2^0^c k'^ 



(A50) 



(A51) 



(A52) 



Gj^„(fc,a;) is not needed in later calculations. The bare 
propagator of 5J±^ is trivially given by 



1 



(A53) 



We finally arrive at the equations of motion as the self- 
consistent form: 



i/i(fc,w) 
rh{k, uj) 
i^(5J||(fc,a;) 



■V;o(fc,a;) 


^SJ°{k,uj)^ 



G"{k,uj)V{k,uj) 



(A54) 
(A55) 

(A56) 
(A57) 

(A58) 
(A59) 



and 

SJ^{k,uj) = dJ°j_{k,Lu) + G'[y±^^(A;,a;), 
where 

4'"{k,uj)^ Gl{k,uj)a^'9r.{k,u), 
6J^^{k,u)^ G|j(A;,a;)6i||(fc,cj), 

dJ°j_(k,u)^ G'i{k,uj)e^{k,u), 

V^^^{k,uj)^-iT,VAk)- [ qXo\q) 

Jqn 

XiP{q)'iJj{k~q). 

Here, ip'^{k,uj), SJ^^{k,uj) and 6J°j_ are the bare variables 
that are the solutions without the nonlinear terms. Iter- 
ating Eqs. (IA54[) and (|A55I) . we can obtain perturbative 
expansions about nonlinear interactions V and V_i_^^. 
We note that the first and third rows of Eq. (jA54|) are 
the equations of motion for the thermal and sound modes, 
respectively, while Eq. (jASSP is for the viscous mode. We 
also stress that Eqs. (|A5ip - (|A53|) are the propagators of 
the thermal, sound and viscous modes, respectively. 
Now, we calculate the two body correlation of ip'^{k, uj) 

and SJj_{k,uj) which are needed in later calculations. 

(^°(fci,wi)^°(fc2,c^2)) =G^(fci,wi)G^,(fe2,^2)ar' 

x{0{ki,uji)0{k2,uj2))- (A60) 



{e{ki,u;,)0{k2, U2)) = 2al\o{k){2nf+H{k^ + fca), 

(A61) 



and 



ujl + \l{k^)xr{ki) 
x(27r)'^+i^(fci +fc2), (A62) 

where 5{ki + k2) = S{ki + k2)S{uji + W2). By a similar 
calculation, we obtain 

{6rAk,,co,)srAk2,c.2))^ 2f1°,?^!-2 (^^(fcl))^.^ 

+ ri^(ki)hc 
x{2nf+^S{ki + k2), (A63) 

where ryo(fe) = Vok^- For a later convenience, we define 

uj^ + mk)xo (k) 
uj^ + T]o{k)hc 



Appendix B: Renormalization of the thermal and 
viscous diffusion modes 



Here, we first derive the RG equations for the thermal 
conductivity and the shear viscosity. Now, we note that 
the sound mode is not a genuine-relevant mode but a sec- 
ondly mode that is strongly affected by order-parameter 
fluctuations but yields only a negligible feedback for the 
order parameters [1^ [sl]. Then, we can neglect the 
sound mode for the minimal critical dynamics; however, 
the bulk viscosity is not renormalized in that case. Here, 
to first analyze the minimal dynamics, we neglect the sec- 
ondly mode, which is renormalized in the next section. 
In that case, the equations of motion are given by 

i^ik.oj) ^^"ik,uj) + Gl{k,io)V^^^{k,ij) (Bl) 



and Eq. (jA55[) . For a diagrammatic treatment , we de- 
note the full and bare variables, the bare propagators and 
the bare correlation functions as Fig. [T] Then, we can 
represent the equations of motion (jBip and (IA55|) as Fig. 

m 

For coarse gaining, we decompose the variables into 
the long- and short-wavelength components as 



'4^ik,uj) = ■0^(fc,a;) -h-0^(fe,w), 



(B2) 



with 



tP^{k,iu)=e{A-SA-k)ip{k,uj), (B3) 
iP^{k,uj) =Q{k - A- SA)'iJj{k,uj), (B4) 

where Q{x) is a step function; i.e., the wavenumber is de- 
composed into < fc < A— (5A and A— 5 A < k < A. Here- 
after, quantities with the suffixes L and S are supposed 
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ip : 



. 



SJ 



SJ° : (?o ■ 



- C"[: -O 



FIG. 1: Diagrams for the full and bare variables, the bare propagators and the bare correlations. 



= \ 



FIG. 2: Diagrams of the equations of motion for the thermal 
and viscous modes. The left and right hand side respectively 
denote Eqs. (|BT|| and (|A55|) . 



to be decomposed as above. To average over the -0°^ and 

~ OS 

6J ^ , wc must solve the equation of motion about them. 
Here, we solve the equations of motion to second order in 

~ ~- OS 

the nonlinear interactions and average over tp^^ and 6J j_ . 
Then, we find the coarse-grained equation of motion for 
■0, which is diagrammatically given by Fig. [3l The last 
two terms in Fig. [3]represent nonlinear interactions being 
of third order, and can be neglected. Furthermore, the 
fifth term vanishes due to the relation between the step 
and delta functions in the loop integral. Introducing the 
self energy S 0^ , which is graphically represented in Fig. 
m we can write the coarse-grained equation of motion for 
tp as 

^L(fc,c.) ^^°^ik,u;) + Gf{k,u;)V^^^ik,u;) 

+0^(fc,w)G^L(fe,w)Sv,^(fe,w). (B5) 
The self energy is given by 



(fc,^) 



{k-V{k-q)-k)xo{q) 



(B6) 



+ Ao(g)xo (q) + Mk - q)h 

where rio{k) = rjok^. Solving Eq. (IBSP about ip^, we 
have 



+ - S^^] V^^j_. (B7) 

where we have used Eq. (jA56l) . Introducing renormalized 
variables as 

{G^R)-\k,u) ^{G''^)-\k,uj) - S^v^(fe,c^), (B8) 
^^^ik,oj) ^G^nik,io)a^^9nik,Lj), (B9) 
we can rewrite Eq. (jA56[) as the renormalized equation of 



motion: 



^i^'^'^{k,co) + GMk,^)V^^^. (BIO) 



We now require that the renormalized propagator has the 
same form as the bare one: 

(G^R)-i(fe,L^) - -iij + \iiD^k''x^\k), (Bll) 

where Ar, is the renormalized thermal conductivity. That 
is, we require that the only transport coefficients are ex- 
plicitly renormalized. The small correction for the ther- 
mal conductivity (5A = Ar — Aq reads 



(5A=- lim [{D^,k^xo{k)) '^T.^^{k,uj)], 

Tc f {k-V{q)-k)xo{q) 
hcD^ Jq Xo{q)xo\q) + Vo{q)hc^ 



(B12) 



We approximate the denominator and the numerator as 

Hq)Xo\q)+Voiq)h-' ^ mik)h-\ (B13) 
Xo\q)=ro + q^ -q^, (B14) 



near the CP [3l|. Then, we find 
Tr f dQ^ 



D^Vo J (27r) 



{k-V{q)-k) 



dqq' 



.d-5 



ASA 



'-^{k-V{q)-k)K'^~HK, (B15) 



D^'no J (27r) 

where d^d is the solid angle in the space dimension d. 
Therefore, we obtain the RG equation for the thermal 
conductivity: 



■A 



Tr 



dn, 



-^{k-V{q)-k)K^-\ (B16) 



9A D^rjiA) J (27r) 

where rjQ is rewritten as »7(A). For the space dimensions, 
d = 4 — e, the angle integral is given by 



—L^k-Viq)-k) = -K,, 



{2n 



(B17) 



where K4 is the surface area of a unit sphere in 4 di- 
mensions divided by {2tt)^. The RG equation in 4 — e 
dimensions reads 



(B18) 
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L S s 

+2X +^^^^^^ 




FIG. 3; Diagrams for the coarse-grained equation of motion for -i/j. Tiie letters, L and S, respectively denote the long- and 
short-wavelength components (see the text below Eq. (|B4|l '). 




where we have introduced 



/(A) 



(B19) 



for a later convenience. 

By making coarse graining of the viscous mode with a 
similar procedures as above, we obtain a small correction 
for the shear viscosity: 

Sv^^ lim [(fc2ft-i(d_i))-i V(Si±(fc,c.))d, 

(B20) 

where {Y,±±{k,ijj))ij is the self energy for the viscous 
mode and given by 

Jq 



-iuj + Xo{q)xo ^(9) + A(fc - q)xo ^{k ~ q) ' 



(B21) 



which is graphically represented as Fig. 2] . In the space 
dimension d = 4 — e, we find the RG equation for the 
shear viscosity 



where the prefactor 1/24 comes from the angular integral 
in Eq. ((B?T1) and the factor (d - 1)-Mn Eq. (|B20)) . 

Differentiating Eq. (|B19p about A, we have the RG 
equation for /(A) 



A^ = (e-i^/(A))/(A). (B23) 



Appendix C: Renormalization of the sound mode 

Next, let us make a coarse graining of the sound mode 
for the renormalized bulk viscosity. Because a feedback 
from the sound mode is neglected, we must renormalize 
the mode with a method separating relevant and sec- 
ondly modes ISSll. Here, we take the method developed 
by Onuki [20l. l26l| , in which RG equations are derived 
from fluctuation-dissipation relations. 

Now, we consider the equation of motion for the sound 
mode, 1X22)) : 



+^2 hcl^ / i^i^Mk - q)] 

Jqfl 



(CI) 
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where the noise term 6'j| satisfies the fluctuation dissipa- 
tion relation: 

x{2T:Y+^5{ki+k2). (C2) 

Since 5J\\ is a conserved density projected onto k, we can 
rewrite Eq. (|C1|) as 

-ioj8J\\{k,uj)^ik-Ii{k,uj)-k, (C3) 

where Ily is the stress tensor. If we take z direction as 
fc, ITzz reads 

+ifczy^/i;:Mj||(A;,w) +^i-^Jfc,cj), (C4) 

where 7r,^j(fc,aj) is the random-stress tensor coming 

from microscopic process and satisfies the relation, ik ■ 

n°{k,iu)-k^ e'^^{k,uj). 

We now consider how Eq. (jCip is affected by the 

coarse-graining procedure. In the coarse-graining proce- 

~ ~ s 

dure, the variables, tp^ , and 6 J are eliminated from 

Eq. (|Cip . The eliminated variables do not disappear 
from the equation of motion but are implicitly contained 
in the noise term. In other words, we convert the macro- 
scopic process in the wavenumber shell A ~ SA < k < A 
into the microscopic process. In this procedure, the noise 
term is implicitly renormalized as follows 

0^{k,Lo) ^ ^{k,uj) + 0^^"°{k,uj), (C5) 

where 

ef''"°{k, io)= ik ■ ^'^'^"°(fc, Lo) ■ k, (C6) 

+l3^'h,^o I i^\q)i^\k - q)] 
+ikvl,h~^5Jl, (C7) 
— T,/32-i/ie7o / ^"(g)^S(fc-g), (C8) 

JqU 

where we neglect the linear terms in Eq. (jC7p that is 
irrelevant for the following argument. The new term 
9^'^"°{k,ijj), being due to the coarse graining, con- 
tributes the transport coefficient through the fluctuation- 
dissipation relation: 

= 2Tckl5v\ki,u:i){2TiY+^5{ki + fca), (C9) 



where we have assumed that the renormalized equation 
of motion has the same form as Eq. (|C1I) . We note that 
this assumption is equivalent to the requirement below 
Eq. (|B10p . Now, we calculate the left-hand side in Eq. 
([^: 

X / (v?(<zi)v;^(fci-gi)vi'(g2)Vi'(fc2-<z2)). (CIO) 

^giOig2Sl2 

Approximating the variable by the bare one, ip^ ~ i^^^ , 
we find 

(0~^'^"°(fcl,L.l)0~^-"°(fe2,C^2)) = (2^)^+M(fci + fc2) 

x2fc2(T,M2-')'7o / Cf{q)Cf{k,-q), (CU) 

where we have used Eq. (jA62p and neglected a term 
corresponding to a disconnected diagram. Then, com- 
paring with Eq. (|C9p . we obtain the correction to the 
longitudinal-kinetic viscosity: 

5v\k,u^) = T^P^^hl^l I Cf{q)Cf{k - q). (C12) 

JqU 

We are not interested in the frequency- or wavenumber- 
dependent bulk viscosity and then take the limit fc, — 
0: 

5z/'= lim 5v\k,uj) 

= T,P2^hl^l I {Cf{q)f. (C13) 

After the integration, we find the RG equation for longi- 
tudinal kinetic viscosity: 

where we have rewritten the static parameter 70 as 7(A) 
to denote its cutoff dependence as mentioned in the text. 
The asymptotic behavior obtained from this RG equa- 
tion is different from the shear viscosity's behavior, so 
we replace above RG equation as 

Although, by this method, we could more easily obtain 
the RG equations for the thermal conductivity and shear 
viscosity, we have taken the diagrammatic method for an 
instructive purpose. 



16 



[1] M. Asakawa and K. Yazaki, Nucl. Phys. A 504, 668 
(1989); 

A. Barducci, R. Casalbuoni, G. Pettini and R. Gatto, 
Phys. Lett. B 231, 463 (1989); Phys. Rev. D 49, 426 
(1994); 

As review articles, see, M. A. Stephanov, Prog. Theor. 
Phys. Suppl. 153, 139 (2004) [Int. J. Mod. Phys. A 20, 
4387 (2005)]. 

[2] M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto, 
Prog. Theor. Phys. 108, 929 (2002). 
Z. Zhang, K. Fukushima and T. Kunihiro, Phys. Rev. D 
79, 014004 (2009). 

Z. Zhang and T. Kunihiro, Phys. Rev. D 80, 014015 
(2009). 

[3] N. Yamamoto, M. Tachibana, T. Hatsuda and G. Baym, 
Phys. Rev. Lett. 97, 122001 (2006); Phys. Rev. D 76, 
074001 (2007). 

[4] M. G. Alford, J. A. Bowers and K. Rajagopal, Phys. Rev. 

D 63, 074016 (2001) 
[5] E. Nakano and T. Tatsumi, Phys. Rev. D 71, 114006 

(2005). 

S. Maedan, Prog. Theor. Phys. 123, 285 (2010). 
M. Buballa and D. Nickel, arXiv:0911.2333 [hep-ph]; 
S. Carignano, D. Nickel and M. Buballa, Phys. Rev. D 
82, 054009 (2010). 

T. Kojo, Y. Hidaka, L. McLerran and R. D. Pisarski, 

Nucl. Phys. A 843, 37 (2010). 
[6] P. de Forcrand and O. Philipsen, JHEP 0701, 077 (2007). 
[7] K. Fukushima, Phys. Rev. D 78, 114019 (2008). 
[8] M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. 

Rev. Lett. 81, 4816 (1998); Phys. Rev. D 60, 114028 

(1999). 

[9] T. Kunihiro, Phys. Lett. B 271, 395 (1991). 

[10] T. Kunihiro, p. 287 in " CONFINEMENT 2000, Proceed- 
ings of International Symposium On Quantum Chromo- 
dynamics (QCD) And Color Confinement (Confinement 
2000), 7-10 Mar 2000, Osaka, Japan,", Edited by H. Sug- 
anuma, M. Fukushima, H. Toki. ( Singapore, World Sci- 
entific, 2001); arXiv:hep-ph/0007173 

[11] Y. Hatta and T. Ikeda, Phys. Rev. D 67, 014028 (2003). 

[12] H.Fujii, Phys. Rev .D 67, 094018 (2003); H.Fujii and 
M.Ohtani, Phys. Rev. D 70, 014016 (2004); H. Fujii and 
M. Ohtani, Prog. Theor. Phys. Suppl. 153, 157 (2004); 
H. Fuju and N. Tanji, J. Phys. G 35, 104060 (2008). 

[13] D.T.Son and M.A.Stephanov, Phys. Rev. D 70, 056001 
(2004). 

[14] M. A. Stephanov, Phys. Rev. Lett. 102, 032301 (2009). 
[15] See, for example, H. Haken, Rep. Prog. Phys. 52, 515 
(1989). 

[16] F.Karsch, D.Kharzeev and K.Tuchin,Phys. Lett. B 663, 



217 (2008). 

[17] P. Romatschke and D. T. Son, Phys. Rev. D 80, 065021 

(2009) arXiv:0903.3946 [hep-ph]]. 

[18] G. D. Moore and O. Saremi, JHEP 0809, 015 (2008) . 

[19] C. Sasaki and K. Redlich, Nucl. Phys. A 832 ,62(2010). 

[20] A. Onuki, Phase Transition Dynamics, (Cambridge Uni- 
versity Press, 2007). 

[21] G.F. Mazenko, NoneqmUbrium Statistical Mechanics, 
(WiLEY-VCH, 2006). 

[22] H. Mori and H.Fujisaka, Prog. Theor. Phys. 49, 764 
(1973). 

[23] Y. Minami and T. Kunihiro, Prog. Theor. Phys. 122, 881 

(2010) . 

[24] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perg- 

amon,New York, 1959). 
[25] B. I. Halperin and P. C. Hohenberg, Phys. Rev. Lett. 19 

(1967), 700; Phys. Rev. 177 (1969), 952; P. C. Hohenberg 

and B. I. Halperin, Rev. Mod. Phys. 49 (1977), 435. 
[26] A. Onuki, Phys. Rev. E 55, 403 (1997) and references 

therein . 

[27] K. Kawasaki, Ann. Phys. 61, 1 (1970). 

[28] K. Kawasaki, Nonequilibnum and Phase transition, 

(Asakura Syoten, 2000), in Japanese. 
[29] H. Mori, Prog. Theor. Phys. 33, 423 (1965). 
[30] B. I. Halperin, P. C. Hohenberg and E. D. Siggia, Phys. 

Rev. B13, 1299 (1976). 
[31] E. D. Siggia, B. I. Halperin and P. C. Hohenberg, Phys. 

Rev. B13, 2110 (1976). 
[32] K. Kawasaki and J. Gunton, Phys. Rev. B 13, 4658 

(1976). 

[33] B.I. Halperin, P.C. Hohenberg, and S. Ma, Phys. Rev. B 

10, 139 (1974); 13, 4119 (1976). 
[34] Y. Akamatsu, T. Hatsuda and T. Hirano, Phys. Rev. C 

79 (2009) 054907 arXiv:0809.1499 [hep-ph]]. 
[35] W. Israel, Ann.Phys.(N.Y.)100, 310 (1976); 

W. Israel, J.M.Stewart, Ann.Phys.(N.Y.)118, 341 

(1979) . 

[36] W. A. Hiscock and L. Lindblom, Phys. Rev. D 31 (1985), 
725. 

[37] C. Nonaka and M. Asakawa, Phys. Rev. C 71, 044904 
(2005) arXiv:nucLth/0410078;. C. E. Aguiar, T. Ko- 
dama, T. Koide and Y. Hama, Braz. J. Phys. 37, 95 
(2007). 

[38] D.M. KroU and J.M. Ruhland, Phys. Lett. 80A, 45 

(1980) . 

[39] M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto, 
Phys. Rev. D 65, 091504 (2002); Prog. Theor. Phys. 114, 
117 (2005). 



